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Q Gravitational radiation reaction affects the dynamics of gravitationally bound bi- 

nary systems via "tail" terms which, at the lowest level, modify the conservative 

m 

dynamics at fourth post-Newtonian order, as it was first computed by Blanchet and 
Q Damour. Here we re-produce this result using effective field theory techniques in 

^ the framework of the closed-time-path formalism. This tail term is the lowest order 

I— I example of a short-distance singularity showing up in the conservative dynamics, 

and it is correctly taken into account within the effective field theory formalism. 
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I. INTRODUCTION 

The forces induced on a isolated system by reaction to the emission of gravitational waves 
and their impact on the motion of gravitationally bound binary systems has been studied in 
great accuracy since the first derivation of the radiation reaction force in General Relativity 
by Burke and Thorne [TH3]. Their phenomenological impact is linked to the forthcoming 
observation runs of the Laser Interferometer Gravitational Observer (LIGO) and Virgo, see 
[1] for the result of a recent search of coalescing binaries, and have already been observed 
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to be at play in binary pulsar systems [5l |6]. 

The motion of coalescing binaries is imprinted in the shape of the emitted gravitational 
waves and the output of gravitational detectors is particularly sensitive to the time varying 
phase of the radiated wave, which has to be determined with high accuracy in order to ensure 
high efficiency of the detection algorithms and faithful source parameter reconstruction. 
The standard approach to describe the motion of coalescing binaries lies within the post- 
Newtonian (PN) approximation to General Relativity, describing the binary system dynam- 
ICS cLS cL perturbative series in terms of the relative velocity of the binary constituents, see 
e.g. [7] for a review. 

The leading effect of radiation reaction modifies the binary dynamics giving rise to a term 
odd under time reversal which changes the energy of the system at 2.5 PN order: this is the 
lowest order at which linear effects of the gravitational radiation enter. 
In [HI [9] the leading non-linear radiation reaction effect has been derived, and it is shown 
to modify the binary dynamics at 4PN order (i.e. at 1.5PN order relatively to the leading 
effect): it belongs to the species of terms dubbed hereditary, as it depends on the entire 
history of the source. In particular it originates from radiation emitted and then scattered 
back into the system by the background curvature generated by the mass M of the binary 
system, hence the name of tail term. Such non-linear 4PN tail term has been shown to 
include both a term non-invariant and a term invariant under time reversal, the latter 
can be incorporated in the conservative dynamics of a binary system, in good agreement 
with computations performed within the framework of the gravitational self-force analysis 
of circular orbits in Schwarzschild background, as found in [TO] . 

Here we present the re-computation of the 4PN tail contribution to the conservative 
dynamics of inspiraling binaries via the use of the Effective Field Theory (EFT) methods 
for gravity introduced in [11]. EFT methods turns out to be useful in problems admitting 
a clear scale separation: in the binary system case we have the size of the compact objects 

the orbital separation r at which the system consists of point particles interacting via 
instantaneous potential, and the gravitational wave-length A (with hierarchy < r ~ 
Ts/v"^ < A ~ r/f, being v the relative velocity between the two bodies) at which the binary 
system behaves as a particle of negligible size endowed with multipoles. 

Using this approach several different groups have re-produced results in the PN analysis 
which have been previously computed in the standard approach both in the conservative 
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Figure 1: Diagram describing the gravitational radiation reaction force at leading order (left) and 
leading non-linear order (right). The thick line represents the massive binary system, the curly 
line the the gravitons emitted and absorbed by the system, the dashed line the potential graviton 
responsible for the Newtonian potential. 

[121 [13] and in the dissipative [H] sector. Moreover new results on the PN analysis have 
been made available by the use of the EFT method in both sectors [IHl [16], as well as in the 
extreme mass ratio limit approach to binary coalescence [17]. We have recently shown [13] 
how to recover the Hamiltonian dynamics at 3PN order using effective field theory methods, 
by implementing an algorithm using the Mathematica [18] software to integrate out short 
range Newtonian gravitons to determine the conservative binary dynamics, paving the way 
to higher order computation not involving tail terms. However at 4PN order, the radiation 
reaction comes into play and we show in the present paper how EFT methods are able to 
recover the 4PN conservative tail, which involves emission and absorption of radiation by 
the gravitational system. 

The leading order radiation reaction effect via effective field theory methods have been 
computed in [19J adopting the closed-time-path formalism, which is the suitable framework 
in which computing amplitudes with time-asymmetric boundary conditions, as the physical 
set up represented by coalescing binaries requires no incoming radiation from the past-null 
infinity. 

In the following we use c = 1 units and the mostly plus signature convention. Contraction 
of space indices are taken with the Kronecker delta. 



II. RADIATION REACTION FROM TAIL TERM 
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Following the EFT framework for non relativistic General Relativity, after integrating out 
the potential gravitons we are left with an effective action at the orbital scale r, describing 
radiation gravitons coupled via moments of the two-particle distribution to a composite 
static object of total mass M given by the sum of the binary constituents masses mi^2- We 
adopt the decomposition of the metric suggested in [22] 

\^ Ai/mpi e-^*/™^' {5ij + aij/mpi) - AiAj/m% J 

The coupling of the radiation gravitons can be read from eq. (24) of [H], here we report the 
relevant terms for the radiation reaction computation, that are 

1 1 



S^ff ^-—Jdr (^M0 + -hAj) , (2) 

with Gn = l/(32vrmp;), R^ioj denotes the appropriate component of the Riemann tensor 
and the moment lij is defined as lij = m„x„jx„j.[26] 

At lowest order, by integrating out the radiation graviton, i.e. by computing the diagram 
in the left of fig. [T| the Burke-Thorne potential is obtained, which has been derived in the 
EFT framework in [T9,j. Corrections to the leading effect appears at relative IPN order 
due to the inclusion of higher multipoles and the IPN modified dynamics of the quadrupole 
[20| |2T] . The genuinely non-linear effect appear at relative 1.5PN order and it is due to 
the diagram in the right of fig. [T| whose calculation is the subject of the remaining of this 
paper. The interactions involved in the 3-particle vertex can be read from eq. (6) of [TH] . 
The effective action obtained by computing the right diagram in fig. [T] is given by the sum of 
six contributions, one for each interaction vertex which can be built out of the metric fields 
0, A and cr, given that the M insertion sources a 0. 

In order to perform the computation of the diagram, boundary conditions asymmetric 
in time have to be imposed, as no incoming radiation at past null infinity is required. 
Technically this is implemented by adopting the closed-time-path (CTP), or in-in formalism 
(first proposed in [23], see [21] for a review) as described in [19]. The CTP construction 
requires to double the field variables, e.g. for a scalar field the generating functional W 
for connected correlation functions has the path integral representation 



(3) 
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It turns out to be useful to perform a change of variables from the 1, 2 base to a +, — base 
defined by 

^- = ^1 - ^2 , ^+ = ^ , (4) 
so that for generic fields "^a, and the doubling of the degrees of freedom gives 
^ ^_ , 

^A^^B ^Al^Bl - ^A2^i?2 = ^A+^B- + , (5) 

For free scalar fields, the path integral in eq. (|3| can be performed exactly, leading to the 
effective action Vr[J+, J_] = i/2 J df^x d^yJD{x)G^^ {x — y)JE{y) which can be written in 
the so-called Keldysh representation [25] (indices D and E taking values + and — ) where 



G^^(t,x) 



^ -zGA(t,x) 

where G~^^ = and Ga{r) is the usual advanced (retarded) propagator 

.S{t + r) 



(6) 



zGAit,^) = ^^(-t)(A+(t,x)-A_(t,x)) = 
-zG'«(t,x) = ^^(t)(A+(t,x)-A_(t,x)) = z 



Airr 

Gnit-t',^) = A+(t,x) + A^(t,x), 



where 



A±(t,x) = j 



-ikx 



(27r)3 2A; ' 

The computation of the right diagram in fig. [l| using eqs. ^ and (|5]), boils down to the 
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following contribution S^'^^^'' to the effective action Seff 



2 V2A J \ A J \ A J J " 47r 



+2 {2a+ardb(l)+dp(f)- + a^ardb(j)+dp(f)+) {6abSrp + SapSrb - Sar^bp) 
+8 daA^r(p+dbA_p {6abSrp — ^ap^rb + ^ar^bp) 
-8 [A+ada(l)+(i)^ + v4_„(9a0+</; 



-320H 



(ti,x)x 



where all possible contractions between fields at different space-time points should be taken 
(ensuring that one of the two fields involved in contraction is at (tj, x)) and are understood, 
an over-dot denotes a time derivative and da = d/dx"", r = |x|, and where 

2 

-^—ij ^ ^ {^■^—ai^+aj ~l~ ^+ai^—aj) ; 

2 



■^=1 (9) 



a=l 



After taking the Wick contractions as summarized in eq. (Al), we obtain mixed terms of 
the type /_/+ and terms quadratic in /_. In order to get the equations of motion for the Xj's 
we have to take the variation with respect to x_j and then evaluate the result at x_j = 0, 
x+j = Xj, so we can discard the terms involving Dh as they are quadratic in J_, as well as 
the terms quadratic in x_ present in in /+. Using the Einstein Hilbert action, the metric 
parametrization in eq. ([T]) and the gauge fixing term as in [13j , the retarded propagators for 
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the 0, A and a fields are 

1 



P[0(t,x)0(t',x)] 



1 



P[A,(t,x)A,(t',x')] = -5. 



> X Gij(t-t',x-x'), (10) 



PK(t,x)a.Kt',x')] = - li5.,S,, + 5u5,,-25.M 



where all other mixed propagators vanish. Introducing the quadrupole moments 

Q±ij = I±ij ~ 7^SijI±kk , (11) 

after performing the calculations reported in the appendix, the final result turns out to be 
MQM) _ ,4^2 f^^r^ u\ t d^Q+iji't') 1 



^^e?r^ = -'^GnM J dtQ..,{t) J ^ dt' ^^e^ . (12) 

This results exhibits a short-distance singularity, which is inherited by the contribution to 
the equations of motion containing the effect of the tail term A^^'^^^Xai- 



A(««*^)x„.(t) = - ^ 
m„ ox- 



(QQM) 



(13) 



The divergence can be regularized by standard methods, i.e. by renormalizing the bare 
source moment Qij to absorb the divergence which then disappears from the physical results 
expressed in terms of the renormalized moment. After integrating by part and regularizing, 
expressing the results in terms of the regularized moment, which we still denote by Qij, one 
obtains: 

A(««*^)x,,(t) = -^x.,(t) dt' ^^^^ log [{t - t')/i] . (14) 

J~oo dt 

where the renormalization scale /i has been introduced. 

The explicit /i-dependence in eq. (14) is absorbed by a compensating dependence on /i in 
the renormalized quadrupole moment. Consistently with renormalization group invariance, 
physical results are independent on the arbitrary renormalization scale /i. In practice one 
can take fi to be of the order of the typical graviton frequency = f/r, or of the typical 
short distance scale in the problem fi = 1/r. Following [10] we adopt fi = 1/r and using the 
trivial identity 
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valid for any arbitrary A, we obtain a piece in the equations of motion which contains no 
exphcit dependence on time and which can be interpreted as a quadrupole-sourced contri- 



bution affecting the conservative motion of the sources Aions'^^^Xai- Using eqs. (14) and (15) 
one obtains 

A(?r ^^a. = l^'^AM^ log Q , (16) 
which coincides with the equations of motion derived in following the result of [8]. 



III. CONCLUSION 



The conservative dynamics of gravitationally bound binary systems is known in literature 
up to the third post-Newtonian order. Such results has been derived by both traditional 
methods and within the context of effective field theory. In both methods long distance 
singularities are met and consistently dealt with. At fourth post-Newtonian order a radiation 
reaction tail term affects the conservative dynamics, involving the regularization of a short 
distance singularity. We have shown that the effective field theory approach consistently 
reproduces the result obtained by traditional methods, adding yet another ingredient to 
support the equivalence of the two approaches and to cross-check old results. 
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Appendix A: Useful formulae 



We give an explicit derivation of the Wick contractions taking from eq. (Isj) to the result 



(14). Within the CTP formalism, in the Keldysh representation, the Wick contractions give 



I [-iGnit - U)\Gh{U - t') + iDuit - U){-iGa{U - 1')) 

+\GH{t - U){-iGa{U - t')) - iGnit - U)h3H{U - t') I-{t') (Al) 

\2{-tGR{t - u)){-tGR{U - m I+{t') 
+/+(t) [2{-iGA{t - U)){-iGa{U - m I^it')]i-tGRiU - 1")) , 



where we have suppressed tensor indices and we have assumed the field R to have non- 
vanishing y4-independent contraction with any of the fields "^a, 4>- In our case we have to 
use the Riemann tensor expanded to first order as 



iOj 



1 1 • 1 • 



(A2) 



where h denotes generically the field 0, A or a. After taking all the Wick contractions in 
eq. ([8| and keeping only the terms linear in /_ one obtains: 



iS 



(QQM) 



-(- 

2 V 2A 



-iM\ /z47r\ 



A 



A J \A{8tt 



dtdt' I dti I drr"^ — | 



I.i,{t)I+ki{t')x 



+325,, ( d, 



l \ 1 



d'',5{t-U-r)dt3t,dk 



5{U -t'-r) 



-8 (^d^^ ^d^t -U- r)d, {d,d, + 5,,dl) {^-^^^-^—^ 

+lQ^dtdida — ^ — dt'dkdb — ^ {SajSu - 5ai5bj + 5ah5ji) 

[0.1) 3,3, (^<i^) 4 (3.3. . «?) i^-^t^ 
--ot, [Oidj + 6ijdt ) I 1 dt^ [dkdi + 5kid^,) ( 



(A3) 
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Using the following 



Air 

J dQ^rirjVkri 



1 ^ij ' 



da 



dadb 



S{t - r) 
r 

d{t - r) 



15 



{SijSki + SikSji + Sii6jk) , 



dadbdc 



6{t - r) 



S{t - r) 



■ii. j-an S{t-r) / Tan 

6{t-r)— [Sab - 3^ 



S{t-r) f Tan 
Oab — <J — ^ 



(A4) 



+3 



8{t - r) 



Sab \-Sac \-Sbc 5 



r 



+3 



6{t - r) ( r, 



n 



Ta r- r-anr-c 



>ab- 



+ Sac \- Sbc 5 



we obtain the following contributions from respectively the aAcj), a4>^, A^(f), A(jP and 
d)^ vertices (we substitute u = t — t' — 2r ): 



■niQQM) _ , G%M 



/oo poo poo -j^ 

dt / dt' I-ij{t)I+ki{t') / dr^ 
-oo J-oo Jo 



X 



5(4) (m) 



32 

T 

16 
15 

16 

y 



6.(^)(.) + 15^ + 18^ + 9^ 



(A5) 



+ 



16 
15 



+ 6 



+ 3- 



u 



5(4) ( 



U 



+ 



15 



+ 5ik5 



ik(Jjl 



SijSki I I 6 



5(5) («) ,J^^\u) ,„5(=')H J^^Ku) 



+ 15 



+ 18- 



+ 9- 
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Now we use 



J, ml 2,.) 



- /6 



(A6) 



to find that the terms contained in the second, third and fifth square brackets of the 
above eq. (A5) (corresponding respectively to the a(f)A, cfxpa and (j)(f)A bulk interactions) are 
identically zero, and in the other lines only the 5^^\u) terms contribute, summing up to 



. n(QQM) 



iG\M- ( 5ik5ji + 5ii5jk - ■7;^ij^. 



ilUjk 



HjUkl 



oo J —c 



dt 



df 



(t - t') 



n(Ar) 



The amplitude index structure projects the symmetric moments onto their traceless part, 
allowing us to replace I±ij with the traceless quadrupole moments Q±ij as expected. 
Note that when integrating by parts the derivatives of the delta's we repeatedly use 

t'=t 



1 



Qit') 



{t - t'Y a - 1 (t - vy-^ 



t'=-oo 



1 



dv [t-vy-^ 



(A8) 



and imposed 



Qit' 



t'=t 







(A9) 



t'=-oo 



a - 1 (t - t'Y-^ 

for any a ^ 1. In dimensional regularization, only the t' = t divergence for a = 1 (i.e. the 
logarithmic divergence) is physical and it introduces the standard l/{d — 3) divergent term 
in the renormalization constant, being d the number of space dimensions. 
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